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ABSTRACT 

Using 2D numerical hydrodynamical simulations of type la supernova remnants (SNR la) we 
show that iron clumps few times denser than the rest of the SN ejecta might form protrusions 
in an otherwise spherical SNR. Such protrusions exist in some SNR la, e.g., SNR 1885 and 
Tycho. Iron clumps are expected to form in the deflagration to detonation explosion model. 
In SNR la where there are two opposite protrusions, termed ‘ears, such as Kepler’s SNR and 
SNR G1.9+0.3, our scenario implies that the dense clumps, or iron bullets, were formed along 
an axis. Such a preferred axis can result from a rotating white dwarf progenitor. If our claim 
holds, this offers an important clue to the SN la explosion scenario. 

Key words: ISM: supernova remnants — supernovae: individual: SNR 1885 — supernovae: 
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1 INTRODUCTION 

Type la Supernovae (SNe la) are thermonuclear detonations of 
carbon-oxygen white dwarfs fWDs: iHovle & Fowled 1 1960b . Six 
SN la scenarios are currently considered, with no consensus on 
even the leading sce nario for SN la. (a ) The core-degenerate 
(CD) scena r io (e.g., Livio & Riessl 120031 : Kash^&Sokeij 1201 ll : 
ISokerll201ll : lllkov & Sokeii 2012112013 : ISoker et aL l2013h . Here 
the WD merges with a hot core of a massive asymptotic gi¬ 
ant branch (AGB) star. Explosion might occur shortly or a long 
time after merger, (b) The double deg enerate (DD) scenario (e.g., 
IWebbinkll 19841 : Ilben & Tutukovll 1 984h . According to this scenario 
a merger of two WDs takes place, but there is no specification of 
the l ater evolution, e.g., how long after merger explosion occurs, 
(e.g.. IPakmor et al.l 2013: Ivan Kerkwiik et aklboiol : iLevanon et alJ 


l2015l : iTanikawa et _al] |20j5|) (c) The 'double-detonation ’ (DDet) 

mechanism (e.g. , IWooslev & Weaver] 1 19941 : iLivne & Arnettl[l995l : 
IShen et al.ll2013h . in which a sub-Chandrasekhar mass WD accu¬ 
mulates a layer of helium-rich material on its surface. The he¬ 
lium layer detonates and leads to a second detonation near the 
center of the CO WD. A variation of this scenario is the triple 
detonation scenario (Papishetal. 20 l j). (d) The single degen¬ 
erate (SD) scenario (e.g., Whelan & Ibenl 1 19731 : iNomotol Il982l : 


[Han & Podsiadlowskil2004r) . In this scenario the WD accretes mass 
from a non-degenerate stellar companion and explodes more or 
less when it reaches the Chandrasek har mass limit, (e) The WD- 
WD collision (WWC) scenario (e.g. iKushnir et al.ll2013l) . In this 
scenario, two WDs collide and immediately ignite, (f) The singly- 
evolved star (SES) scenario (IChiosi et alJl2015) . In this scenario it 
is claimed that pycnonuclear reaction might be able to drive pow¬ 
erful detonations in single CO white dwarfs. 

A comparison between the first fi ve vari o us sce n arios is given 
in ITsebrenko & Sokej d2015ah and ISokerl d2015h . IPapish et all 


( 2015 ) list two problems with the SES scenario. When confronting 
the scenarios with observations, as well as when trying to account 
for a specific observational property, it is imperative to consider all 
scenarios. 

In the present study we focus on the observational feature of 
two opposite ’ears’ in some SN la remnants (SNRs la). An ‘ear’ is 
defined as a protrusion fro m an otherwise round SNR surface. 

In a previous paper (iTsebrenko & Sokei] 120131) we showed 
that opposite ears observed in some SNRs la might be attributed 
to two opposite jets working either a long time before the explo¬ 
sion or very shortly (within hours) before the SN explosion. In 
the first possibility the ears are formed in the circumstellar matter 
(CSM) that was ejected a long time before explosion, much as ears 
formed in planetary nebulae (PNe). The genera l structure of a SN 
la inside a PN was mentioned bv lDickel & Jonesl (Il985l) . although 
they did not refer to ears or to je ts. The formation of ears in this 
case of SN la inside PNe (SNIP; ITsebrenko & Soked 12015 ah was 
simulated before for the Kepler SNR (iTsebrenko & Sokerl 120 1 3h 
and SNR G 1.9+0.3 (ITsebrenko & Sokerl l2015bl) . This type of in¬ 
teraction can occur in the SD or CD scenarios. But if the CSM 
i s massive, > 0.1 M p, then only the CD scenario can work. In 
ITsebrenko & Soke J ( 2013) we also simulated jets launched by an 
accretion disk around the exploding WD very shortly before the ex¬ 
plosion. Such disks ca n be formed from the merger of two WDs in 
the DD scenario (e.g., Yoon et al1l2007l : iLoren-Aguilar et al1l2009l : 


I.Ti et al.l201 3l : Izhu et al.fcoN) . or in the CD scenario. Another pos¬ 

sibility to for m a small ear in the SNR is by the instabilities arising 
in the ejecta (iWarren & Blondinll2013l) . 

In the current paper we continue our study of ears in SNRs, 
but study their formation by dense clumps, or ’bullets’, along 
an axis inside the ejecta from the SN. We would like to iso¬ 
late the effects of iron bullets formed in the explosion from pos¬ 
sible effects of dense CSM. We therefore differ from our previ- 
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ous works on ears formation result ing from ejecta-CSM interaction 
(iTsebrenko & SokeJl2013Ll2015blh and study an SN explosion in¬ 
side a low-density ISM/CSM. 

We are motivated by the observations of iron clu mps a long 
two almost perpendicular axes in the SNR 1885 dFesen et al.l2014l) . 
and by the pres ence of heavy metals inside the ears in Tycho 
JChiotellis et alJEoTs ). We note though that in Tycho’s SNR the 
ears are not two opposite ears. Such nickel/iron clumps might 
form during the deflagr ation phase of the defl a gration to deto¬ 
nation (DTD) scenario JSeitenzahl et al.112013h . lSeitenzahl et al.1 
(120131) performed 3D simulations of the delayed-detonation explo¬ 
sion model, and found the stable neutron-rich iron-group isotopes 
at intermediate velocities in a shell surrounding a nickel-rich core. 
We are mot ivated by thi s finding as well, but we note that the recent 
analysis bv lMazzali et al.1 d2015h of the first year spectral evolution 
of the SN 201 lfe does not support the claim for neutron-rich iron 
group isotopes moving outside the nickel c ore. SN e ject a clum ps in 
core collapse SNRs have been simulated bv lMiceli et al.1 ( 20131) . ex¬ 
plaining observed protrusion features in the Vela SNR. The role of 
ejecta clumping in forming protrusions in SNRs la has been stud¬ 
ied also bv lOrlando et al.1 l2012h for SNR 1006. We perform 2D 
hydrodynamical simulations of an SN la explosion already hav¬ 
ing an iron clump or a ’bullet’ inside the SN ejecta, and follow the 
evolution of the SNR. 

The 2D numerical code and initial setting and conditions are 
described in Section [2 In Section [3] we show that the ’ears’ in the 
otherwise round shapes of SNRs la can be explained by the pres¬ 
ence of dense iron clumps in the outer ejecta. Our discussion and a 
short summary are given in Section [4] 


2 NUMERICAL SETUP 

The simulations are performed using the high-r esolu tion multidi¬ 
mensional hydrodynamics code FLASH 4.2.2 dFrvxell et alJl2000h . 
We employ a 2D adaptive mesh refinement grid having axisym- 
metric cylindrical (x, y ) geometry, with the y -axis serving as the 
rotational symmetry axis. The length of each axis is A = 20 pc. 
The SN ejecta density p e jecta is modelled by an exponential density 
profile (iDwarkadas & Chevalier! 19981) . 

Pejecta — Aexp( u/u e jecta)^ •> (1) 

where u e j e cta is a constant which depends on the mass and kinetic 
energy of the ejecta, 

Ejecta = 2.44 x lO 8 ^ 2 ' cm s -1 , (2) 

Mch H 1.4M 0 , F? 5 i is the explosion energy in units of 10 51 erg, 
and A is a parameter given by 

A = 7.67 x 10 6 ' E ~ 1 3/2 g s 3 cm" 3 . (3) 

We take the parameters as in ITsebrenko & SokeJ (120131 . 12015bh : 
exploding mass of Mch = 1.4M©, explosion energy of 10 51 erg, 
and maximum ejecta velocity of usNm = 20, 000 km s -1 . 

Our simulations start around 16 yrs after the SN explosion. 
By this time, the SN ejecta has reached a distance of 0.32 pc from 
the explosion location. The mass of the iron clump is taken to be ~ 
0.01 — 0.02M©, and its density pbiob is taken to be 1.5-3 times the 
density of the adjacent SN ejecta, p e jecta (equation [I}. The length 
of the clump is taken to be 0.13 pc and its width (twice the distance 
from the symmetry axis) is taken to be D = 0.02 — 0.03 pc. We 
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Figure 1. The initial geometry of the supernova ejecta and the iron clump. 
Color represents the density according to the color-bar on the right. The 
initial iron clump density at each radius is taken to be 3 times the density of 
the ejecta at the same radius. Blue thin line marks the surface of the clump. 
The total mass in one clump is ~ O.OIM©. The mass of the ejecta inside a 
similar volume would be a third of that mass. 


note that lOrlando et al. ( 2012 ) populate the SN ejecta with many 
much smaller clumps having size of 0.01 — 0.02 pc and then let the 
ejecta evolve, starting the simulation 10 yrs after the SN explosion. 
We start with a smooth SN ejecta having a single clump on each 
side along the symmetry axis (we simulate only one side). We 
take the initial ISM density outside the front of the ejecta to be 
Pism = 10 -24 g cm -3 . The initial density setup is presented in 
Figure [TJ and the initial location of the iron blob is marked by a 
thin blue line. The initial velocity of the material in the clump is 
identical to the velocity of the adjacent SN ejecta at the same radial 
distance, which expands radially with velocities up to usNm- 


3 RESULTS 

As the SN ejecta interacts with the ISM, part of the kinetic energy 
of the ejecta is deposited in the interstellar medium (ISM) and the 
ISM is accelerated outward, as the SN ejecta is decelerated. The 
denser iron blob continues to advance into the ISM, and forms an 
ear in the shape of the resulting SNR. Two shock waves are formed, 
one running outward into the ISM and one inward into the incom¬ 
ing ejecta. A contact discontinuity develops between the shocked 
ejecta and the shocked ISM. These, as well as various physical pa¬ 
rameters, for our fiducial run are presented in Figure [2] Panel [2a] 
shows the density structure, in which the iron clump material is 
outlined by a thin blue line. The time shown is years after sim¬ 
ulation starts. We present results at t — 1585 yrs when the for¬ 
ward shock has reached a radius of r s h oc k = 6.9 pc (beside the 
iron clump that is further out). The total swept-up ISM mass is 
Mism(£ = 1585) ~ 20M©. The ejecta has been substantially de¬ 
celerated. The ear feature is clearly apparent. 

The relevant flow properties are as follows, (z) Instabilities. 
As the ISM decelerates the ejecta a pressure gradient with decreas¬ 
ing pressure toward the center is developed within the ejecta. The 
density on the other hand increases from the ISM toward the center. 
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This situation is prone to Rayleigh-Taylor (RT) instability. Indeed, 
well developed RT fingers and blobs are clearly seen in the den¬ 
sity and temperature maps. However, these do not change the over 
all shape of the SNR and its boundary with the undisturbed (pre¬ 
shock) ISM. The instabil ities are similar to the results obtained by 
IWarren & Blondinl J2013I) . 

(ii) As expected, the denser iron blob deceleration is lower, 
and it penetrates the ISM. An Ear structure has emerged, hav¬ 
ing mostly material originally in the Fe bullet. This feature re- 
sembl es morpholog ically the FeL emitting Ear structure in Tycho’s 
SNR JChiotellis et al .1120131) . A quantitative comparison between 
this feature in our simulation and the one in the X-ray image of Ty¬ 
chos SNR will require further numerical simulation and deeper ob¬ 
servational analysis. These are beyond the scope of this exploratory 
study. Our result is similar to the protrusion obtained by the sim¬ 
ulation of the Vela SNR performed bv lMiceli et ali ( 20131) . The ear 
feature in our simulation is much large r than the many slight pertur¬ 
bations to SNR sphericity obtained bv I Orlando et all (120121) . (Hi) 
A vortex has developed relative to the rest frame of the moving iron 
clump. Such a vortex can increase the survivability of the clumps 
as it is being decelerated by the ISM, making the Ear a long lasting 
structure. To emphasize the vortex we have subtracted a velocity of 
(v x ,v y ) = (1600, 200) km s -1 from all points. The color coding 
in the velocity panel gives the absolute magnitude of the velocity 
before this subtraction. 

Figure [3] shows our fiducial run at an earlier evolution time, 
t=475 yr. By this time the mass of the swept ISM is comparable 
with the SN ejecta mass. The ear feature is distinctly visible at this 
early evolution time as well. 

To check the robustness of our results, we varied two param¬ 
eters: (z) the iron blob density; (ii) the iron blob width along the 
x axis. Figure |4] shows the flow properties for an iron blob having 
a density of 1.5 times the density of the adjacent SN ejecta, and 
width of D — 0.02 pc as in the fiducial run. Figure \5\ shows the 
flow properties for an iron blob having a density of 3 times the den¬ 
sity of the adjacent SN ejecta, but a wider blob with D = 0.03 pc. 

Figures |4| and [5] show that varying some of the iron blob pa¬ 
rameters does not affect the presence of the ear feature, only its 
size. Overall, our results show that an iron blob moving alongside 
the SN ejecta with similar velocity but slightly higher density (1.5- 
3 times the SN ejecta density) can form a noticeable protrusion 
from the otherwise spherical shape of the SNR. This might be the 
case in SNR 1885, having iron tongues (iFesen et al.ll2014r) and Ty¬ 
cho^ S NR, where m etals are present inside the ears protrusions 
dChiotellis et al.ll2013l) . 


4 SUMMARY 


Although many SNRs la exhibit almost spherical large-scale sym¬ 
metry, there are several SNRs la that show a deviation from spher¬ 
ical large-scale str ucture. Single prot rusions exis t in some SNR 
la, e. g., SNR 1885 jFesen et al.ll2014l) and Tycho dChiotellis et al.1 
l2013h . while two opposite protru s ions, or ’ears’, are present 
in Kepler’s SNR (Reynolds et al. I 120071) and SNR Gl.9+0.3 
dRevnolds et al.ll2008l) . We examined a possible scenario to cre¬ 
ate such protrusions, with 2D numerical simulations of SNR la, 
as described in Section [2] The initial set-up includes a dense bul¬ 
let, supposed to be iron-bullet, moving with the rest of the ejecta as 
described in Fig. Q] Such iron clumps are expected to form in the 
deflagration to detonation explosion model dSeitenzahl et al.l20l3h . 
Our results are presented in section [3] Figs. [2| and [3] show the flow 


parameters for our fiducial case, having a blob density 3 times the 
density of the adjacent SN ejecta and blob width D = 0.02 pc, 
at two times. In Fig. |4|we presented the results for a case with a 
blob having only 1.5 times the density of the adjacent SN ejecta. In 
Fig. [5] we presented the results where the blob density is as in the 
fiducial run, but twice as wide. In all cases a distinct protrusion in 
the SNR shell was obtained. 

In SNRs la having two opposite protrusions, such as Ke¬ 
pler’s SNR and SNR Gl.9+0.3, our scenario implies that the dense 
clumps were formed along a common axis. Such a preferred axis 
can result from a rotating white dwarf progenitor. Exploring the 
mechanism that can produce such axis-directed ’bullets’ ejecta will 
be at the focus of a future work, while the present study is limited 
to checking the effect of such ’bullets’ on the shape of the SNR. At 
this point we merely take the axis-oriented clump ejecta as a work 
hypothesis. If our claim holds, this can offer an important clue to 
the SN la explosion scenario, as it requires the white dwarf pro¬ 
genitor to possess a preferred axis. Namely, it might require a 
rapidly rotating white dwarf progenitor. 

The proposed process forming tw o opposite ears is a third al - 
ternative to our previous two processes (iTsebrenko & Sokerll 20131) . 
One is the launching of jets a long time before explosion into 
a circumstellar matter (CSM) shell, namely, in a planetary neb¬ 
ulae phase. The SN la explosion then is termed SNIR In the 
second alternative the two opposite jets are launched shortly be¬ 
fore explosion. In the presently studied process the ear feature 
is obtained without having a CSM shell. It requires though, that 
the SNR interacts with a substantial mass of ISM, such that the 
ejecta is slowed down. The proposed iron bullets scenario studied 
here could be an alternative process to our previously suggested 
ideas, lo wering the percentage of SN e la that originated inside PNs 
(SNIPs) (ITsebrenko & SokeJl2015al) . But still, a rapidly rotating 
white dwarf progenitor is required. 
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Figure 2. The flow properties at t = 1585 yrs for our fiducial case with a blob width of D = 0.02 pc and blob density of pblob( r ) — 3p e jecta(r), whose 
initial setting is shown in Fig. [I], (a) Density map. The blue contour line shows the location of material initially located within the iron clump. The dashed 
black contour line marks the contact discontinuity, i.e., the boundary between the shocked ejecta and the chocked ISM. (b) Temperature; (c) Pressure; (d) 
Velocity magnitude according to the color scheme. The flow direction and magnitude as given by the black arrows are calculated after subtracting 200 km s -1 
from the x velocity component and 1, 600 km s -1 from the y velocity component. This was done to emphasize the vortex present around (0.2, 6) pc. 
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Figure 3. Same as Figure [2] but at an earlier time (t=475 yr). 
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Figure 4. Same as Figure[2] but the initial density of the iron clump is Pblob(r) = 1 .5p e j e cta(r). 
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Figure 5. Same as Figured but the initial blob width is D = 0.03 pc. 
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